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Motivated by the first experimental evidence of ferromagnetic behavior in a three-dimensional 
ultracold atomic gas, we explore the possibility of itinerant ferromagnetism in a trapped two- 
dimensional atomic gas. Firstly, we develop a formalism that demonstrates how quantum fluc- 
tuations drive the ferromagnetic reconstruction first order, and consider the consequences of an 
imposed population imbalance. Secondly, we adapt this formalism to elucidate the key experimen- 
tal signatures of ferromagnetism in a realistic trapped geometry. 
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I. INTRODUCTION 



Itinerant ferromagnetism is a ubiquitous strongly cor- 
related phase of matter in the solid state. The theo- 
retical study of itinerant ferromagnetism dates back to 
the pioneering work of Stoner [l[ and Wohlfarth, which 
showed that ferromagnetism emerges as repulsive pair- 
wise interactions between electrons overcome the kinetic 
energy penalty of polarization. Subsequent theoretical 
work has determined that soft transverse magnetic fluc- 
tuations have the potential to drive the ferromagnetic 
transition first order before the quantum critical point 
is reached (sHl]- Phenomena consistent with a first or- 
der transition have been observed in the solid state; 
though it is difficult to determine whether they are due 
to soft magnetic fluctuations or the coupling of the mag- 
netic moment to phonon degrees of freedom. However, 
Jo et al. J9| have recently presented the first tentative 
evidence jloL lllj of itinerant ferromagnetism in an ultra- 
cold atomic gas. The cold atom gas is a clean system 
in which to study ferromagnetism, completely devoid of 
the interfering phonon degrees of freedom encountered in 
the solid state, so gifts investigators with a valuable tool 
with which to answer long-standing questions about solid 
state ferromagnets. Furthermore, ultracold atoms exper- 
iments also present a unique opportunity to explore fun- 
damentally new physics associated with ferromagnetism 
including the consequences of population imbalance Q, 
a conserved net magnetization [12[ , the damping of fluc- 
tuations by three-body loss 13[ , spin drag [14j , and mass 
imbalance. Here we aim to take advantage of the high 
levels of control investigators can exercise over the exter- 
nal potential trapping the gas and turn to study ferro- 
magnetism in a two-dimensional thin film. 



Itinerant ferromagnetism is difficult to observe in two 
dimensions in the solid state [HI, However, it could 
be realized in an ultracold atom gas by using counter- 
propagating lasers to create one-dimensional potential 
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which will lead to a stacked two-dimensional gas. The 
system also offers investigators the opportunity to study 
the possibility for a superconducting instability to emerge 
near to the ferromagnetic phase transition [17| ■ The two- 
dimensional system is of particular interest in this case as 
it could shed light on high temperature superconductivity 
where antiferromagnctism competes with d-wave super- 
conductivity to form the ground state. Here we adapt the 
formalism introduced for the three-dimensional case Q to 
expose the contrasting behavior of the two-dimensional 
ferromagnet. We develop a formalism that captures the 
effects of transverse quantum fluctuations and explore 
how they renormalize the effective interaction strength. 
We then address how population imbalance modifies the 
behavior of the atomic gas before studying ferromagnetic 
ordering in a trapped geometry. 



II. FIELD INTEGRAL FORMALISM 

It has been long established that quantum fluctuations 
in a three-dimensional fcrmionic gas with repulsive in- 
teractions have the potential to drive the ferromagnetic 
transition first order (3411 ■ To investigate the impact 
of quantum fluctuations in a two-dimensional fermionic 
gas we explore ferromagnetic reconstruction within the 
setting of an atomic gas, adapting the phenomenology 
developed for the three-dimensional case in Ref. 0- We 
adopt this formalism because unlike the Eliashberg the- 
ory Q it provides an exact expression for the free energy 
which then allows us to make a prediction of the criti- 
cal interaction strength for the onset of ferromagnetism 
and study the atomic gas within a harmonic well. More- 
over, ab initio Quantum Monte Carlo calculations p| [Tc| 
have recently been used to verify the three-dimensional 
formalism, which should therefore provide a solid foun- 
dation from which to study the two-dimensional case. 
Although the atoms do not carry spin, we discriminate 
between the two fermionic species with a pseudospin 
(7 E The species cannot interconvert so separate 

chemical potentials \x a tune the population imbalance, 
which in turn pins the net polarization along the pseu- 
dospin direction. However, when the spontaneous mag- 



2 



nctization formed exceeds the population imbalance, a 
nonzero in-plane magnetization emerges. To study the 
potential for ferromagnetic ordering we express the par- 
tition function as a fermionic coherent state path integral 
Z = Tre-^-^) = / V^c~ s with the action 



action to second order in g leads to 
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J P dr J d 2 r with reduced temperature f3 = 
and e^. denotes the dispersion. As we wish 
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to investigate two-dimensional fcrromagnctism we have 
constrained the spatial integral to a plane. A two- 
dimensional atomic gas could be realized experimentally 
using counter-propagating laser beams whose antinodes 
at half-wavelength spacing b will define stacked quasi- 
two-dimensional layers. Though at finite temperature the 
ferromagnetic ordering is only marginally stable, long- 
range order should be stabilized by the weak inter-plane 
coupling fl2| . The repulsive contact interaction param- 
eter g = g5 3 (r) that can be tuned with a Feshbach res- 
onance [2| is linked to the s-wave scattering length a in 
three dimensions through g = y/2/na/b [23| . Unique to 
two dimensions, the interaction strength is independent 
of density. This means that within a trapped geometry 
the entire atomic gas will experience the same effective 
interaction strength and therefore adopt the same polar- 
ization. 

To develop an effective Landau theory of the mag- 
netic transition, Hertz introduced a scalar Hubbard- 
Stratonovich decoupling of the two-body interaction term 
in the spin channel [2j. However, this form of de- 
coupling neglects the potential impact of soft trans- 
verse field fluctuations, which in three dimensions are 
responsible for driving the second order transition first 
order 0, [!]• Therefore, we will introduce a general 
Hubbard-Stratonovich decoupling that incorporates fluc- 
tuations in all of the spin <f> and charge p sectors. In- 
tegrating over the fcrmion degrees of freedom yields 
Z = J cT s T><$fDp with the action 



,5' 



Tr\n[d T + ec-p^+gp-gcr-4>]. (2) 



At this stage a saddle point analysis would determine the 
mean-field values of p and (p. However, quadratic fluctua- 
tions in these auxiliary fields rcnormalizc these equations. 
Therefore we introduce the putative saddle point values 
po for density and m for magnetization, integrate out 
fluctuations in the auxiliary fields, and finally minimize 
the energy to determine po and m. It is also convenient 
to rotate the z-axis from the quantization axis to lie along 
the direction of the saddle point magnetization m, with 
components labeled by s G { + ,—}. After integrating 
over fluctuations in both the density p and magnetiza- 
tion channels (p to Gaussian order, an expansion of the 
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where we have defined the spin-dependent polarization 



operator n s 
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g\m\. The contact interaction means that an unphysical 
ultraviolet divergence arises from the term in the action 
that is second order in g. To remove it we must affect 
the standard regularization of the linear term g (to 2 — p 2 ) , 



setting g 



£k 2 — e k 3 — £ k 4 ) 1 Oil, where the prime indicates that 
the summation is subject to the momentum conservation 
condition ki +k2 = k3 + k4, and A denotes the total area 
of one stacked layer. 

Finally, after carrying out the remaining Matsubara 
summations, one obtains the following expression for the 
free energy: 
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where n s {e) = !/[! + e f3 ^^^ s \ m \VV^ a / b )] is the Fermi 
distribution, and N s = X)k ns ( £k )- To evaluate the final 
nine-dimensional integral in Eq. Q numerically we em- 
ploy the re-parameterization outlined in App.lAlto reduce 
it to a four-dimensional integral. Moreover, as we are in- 
terested in searching for cxtrcma in the free energy with 
changing polarization we can differentiate our expression 
with respect to magnetization, which at zero temperature 
further reduces the integral to just three dimensions. 

To highlight the potential importance of fluctuation 
corrections we briefly study the contribution to the en- 
ergy from particle-hole excitations around momentum 
2kp . At zero temperature a non-analytic contribution to 
the free energy of the form |to| 3 In to 2 emerges. The same 
non-analyticity was found diagrammatically in Refs. 
The formation of a finite magnetization increases the 
phase-space available for the formation of virtual inter- 
mediate pairs of particle-hole pairs, and this phase space 
enhancement donates a non-analytic term to the free en- 
ergy giving the transition the potential for first order 
character. In the next section we study the effect that 
this non-analyticity has on the phase diagram. 



III. PHASE BEHAVIOR 

With the formal development of the theory complete 
we will now apply the formalism to explore the implica- 
tions of fcrromagnctism in the two-dimensional atomic 
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FIG. 1: (Color online) (a) The growth of magnetization m 
with scattering length for different temperatures. The in- 
set figures show the energy landscape with magnetization for 
T = either side of the first order transition, (b) the phase 
diagram of temperature with scattering length shows the first 
order (dashed red line) and second order (solid red line) 
(quasi)ferromagnetic ordering from the paramagnetic phase. 



FIG. 2: (Color online) (a) The growth of magnetization m 
with scattering length at T — 0. (b) the T = phase diagram 
for imposed population imbalance p with scattering length 
a/b shows the first order (dashed red line) and second order 
(solid red line) ferromagnetic ordering from the unmagnetized 
(UnM) to the partially magnetized (PM) and fully magnetized 
(FM) regions, here unmagnetized refers to having no in-plane 
magnetization. 



Fermi gas, and critically compare the results with the 
three-dimensional case [7[- Before we study the phase 
diagram of the fluctuation corrected free energy, to make 
contact with the conventional Stoner theory we first con- 
sider the result of a direct saddle point approximation 
scheme in which the second order term in the free en- 
ergy is neglected. In this approximation at zero temper- 
ature the free energy is F = (1 + a/b\j2-K i l 2 )ix 2 /2ir + 
(1 — a/b\/2ir 3 / 2 )m 2 /27r/x 2 . This expression is exact, and 
with magnetization featuring only as the lowest available 
term in a Landau expansion its analysis is straightfor- 
ward. For a < \f2-K 3 l 2 b ps 7.8746 this model predicts 
that the gas is paramagnetic, whereas for a > \f2 r K i l 2 b 
the system is fully polarized, a scenario that remains un- 
altered with the introduction of population imbalance. 
An immediate corollary is that the spontaneous magne- 
tization formed is independent of the local density, which 
also holds true when fluctuation corrections are taken 
into account. Therefore, within a trap, the entire atomic 
gas adopts the same polarization. 

Having studied the mean- field limit, we now consider 
the repercussions of fluctuation corrections on the be- 
havior of the magnetization. To orient our discussion, 
we first consider a gas with equal populations of up and 
down-spin atoms. As shown in Fig. [ljai) at scatter- 
ing lengths below a ps 3.9456 the energy profile pos- 
sesses a single minimum at zero magnetization. With 
rising interaction strength a second minimum in the en- 
ergy landscape develops at m ps 0.6, which, with ris- 
ing scattering length, deepens in Fig. [TJaii) to become 
the global minimum at a ss 3.9536 and m ps 0.8. At 



this scattering length the system undergoes a first order 
transition from m — into the polarized regime with 
m ps 0.8. As shown in Fig. QJa) with a further increase 
in the interaction strength the magnetization saturates 
at a scattering length a ps 4.0486. Fluctuation correc- 
tions have had significant impact: they have driven the 
ferromagnetic transition to a significantly weaker inter- 
action strength (a ps 3.9536) compared to the mean-field 
case (a ps 7.8746). At this weaker interaction strength 
the m 2 term in the free energy has a positive coefficient, 
and the ordering is driven by the non-analytic |m| 3 lnm 2 
term. The abetment of the transition by fluctuation cor- 
rections and reduction in interaction strength at which 
fcrromagnctism is seen is common to both the two and 
three-dimensional cases, though in two dimensions the 
transition is immediately to full polarization at mean- 
field level and fluctuation corrections drive a first order 
transition at a weaker interaction strength. 

We now turn to address the behavior of the phase tran- 
sition at finite temperature in Fig.[TJa). Increasing tem- 
perature dulls the fluctuation corrections and the scatter- 
ing length of the first order transition rises and the mag- 
netization following the transition is reduced. Fig. [lib) 
shows that at T ps 0.28Tp a tricritical point emerges 
and the system reverts to second order behavior. The 
Mermin-Wagner-Hohenberg theorem [2(| states that al- 
though an ordered phase can exist in two dimensions at 
zero temperature, at any finite temperature fluctuations 
will destroy long range correlations in the system, and 
the state will be characterized by exponentially decaying 
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correlation functions. Therefore we denote the ferromag- 
netic state as a "quasiferromagnet" (ferromagnet with 
fluctuating polarization direction). However, since the 
two-dimensional gas is experimentally realized in a se- 
ries of disks, each one can couple to its neig hbors and 
tunneling should stabilize the phase [HI, H3 . Further- 
more, the Mcrmin-Wagncr-Hohcnberg theorem is valid 
only in the thermodynamic limit and does not apply to 
finite-sized systems. For a two-dimensional Bosc gas with 
attractive interactions it has been shown that a poten- 
tial trap restricts the system and stabilizes a quasi-Bosc 
Einstein condensate [2l|. In a similar way the harmonic 
trap should stabilize a ferromagnetic phase. So far we 
have focused on how equilibrium properties can stabilize 
the ferromagnetic phase, however, there are also non- 
equilibrium aspects to consider. Within the current ex- 
perimental realization of cold atom gas ferromagnetism 
three-body losses necessitate that the experiment be per- 
formed out of equilibrium. Following a quench small 
ferromagnetic domains are formed [22[ which then grow 
steadily The final size of these ferromagnetic do- 

mains ~ 6/fc F M 

at T = O.lTp and kpa = 2 is small 
compared to the length-scale of the thermal fluctuations 
given by aexp[27r(2fc F a/7r- 1)T F /T] w 10 7 /fc F This 
means that at sufficiently low temperature fluctuations 
will not disrupt the ferromagnetic state and so in exper- 
iments a true ferromagnetic phase should be observed as 
shown in Fig.[Tl 

Having addressed the situation without population im- 
balance, we now consider how a fixed spin population 
imbalance influences the phase diagram. The two con- 
stituent species cannot interconvert so an initial popu- 
lation imbalance is maintained by the difference in their 
chemical potentials. However, if energetically favorable, 
the gas can become more polarized either by phase sep- 
aration or the development of an in-plane magnetic mo- 
ment. As shown in Fig. [IJai and aii), at weak inter- 
actions such that a ;$ 3.9456, the energy monotonically 
increases with magnetization, but the magnetization re- 
mains pinned to the minimum value defined by the pop- 
ulation imbalance p. However, with rising interaction 
strength a second minimum develops in the free energy 
landscape from a w 3.9456 and m sa 0.6. If that magneti- 
zation exceeds the population imbalance, then as shown 
in Fig. &b) the system will phase separate between this 
minimum and that at zero magnetization, with relative 
fractions governed by the Maxwell construction. When 
the emerging minimum becomes the global minimum at 
m pa 0.8, then gases with a lower population imbalance 
enter this global minimum with an appropriate in-plane 
magnetic moment. As the magnetization of the mini- 
mum rises it envelops systems with higher population 
imbalance, and tracks the magnetization curve shown in 
Fig.[5ta) until it reaches full polarization at a/6 w 4.048. 
Like the three-dimensional case, the population imbal- 
ance renders the characteristic interaction strength of 
the transition to be almost constant up to an imposed 
population imbalance of p « 0.8, which could be a key 



experimental signature of first order behavior. 



IV. TRAPPED GEOMETRY 

Having addressed the phase behavior of a uniform sys- 
tem, to make contact with the experiment we now turn 
to address the atomic gas trapped within the spherical 
potential V(r) = uir 2 /2. Following the program devel- 
oped in Refs. [l0l - fl2T ] we aim to minimize the free en- 
ergy within the local density approximation using the 
kernel /(r) = F(r) + V(r)[n + (r) + n_(r)] - 7+n+(r) + 
7_n_(r), here F(r) denotes the energy kernel Eq. dU) 
evaluated with the local chemical potential at r. The 
Lagrange multipliers "f± enforce the constraints of con- 
stant number of atoms imposed by the trap geometry 
Ntot = J[ n +{ r ) + n_(r)]d 2 r and population imbalance 
P < J[ n +( r ) ~ n-(r)]d 2 r/N to t', without loss of generality 
we assume that p > and therefore 7+ > 7_. To study 
the effects of spatial density variations we invoke a local 
density approximation that enables the variational min- 
imization 5f/6n s (r) and yields the simultaneous equa- 
tions for the effective local chemical potentials £t±(r) for 
the species in the rotated spin basis 
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These equations can be understood as having been con- 
structed out of three orders of perturbation theory. The 
lowest, independent of the scattering length a, corre- 
sponds to the Thomas-Fermi approximation within the 
confining potential, the term first order in a introduces 
the mean-field energy penalty of the interaction whereas 
the second order term introduces the energy associated 
with magnetic quantum fluctuations. The detailed study 
of the uniform system revealed that the polarization de- 
pends only on the interaction strength and not spatial 
density variations, meaning that the ratio n+(r)/n_(r) 
and therefore /i+ (r) (r) is constant across the trap. 
Therefore, the two equations reduce to just one that is 
solved by iteration. 



A. Heuristic observations 

We first address what can be determined about the 
behavior of the atomic gas heuristically before present- 
ing the results of the full solution of Eq. ([5]). To de- 
velop our intuition we focus on perhaps the most physical 
quantity that can be measured by experiment, namely 
the cloud size. To start the analysis we consider the 
non-interacting limit a = where the system is un- 
polarizcd and the effective chemical potentials given by 
Eq. follow the familiar Thomas- Fermi form. The root 
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mean square (RMS) radius would increase with popula- 
tion imbalance as [l + lyq^) 3 ^ 2 ] 1 ^ 2 /"^ due to the increas- 
ing Fermi degeneracy pressure. With weak interactions 
a <C v / 2tt 3//2 & we need consider Eq. ([5]) only to first or- 
der in a which yields /i±(r) = 7± — ^( r ) — amax[7zp — 
y(r), 0]/2 1 / 2 7r 3 / 2 6. The first order term reduces the ef- 
fective chemical potential so to conserve the total number 
of trapped atoms we renormalize the Lagrange multipli- 
ers upwards from the Thomas Fermi value by a factor of 
1 + a(l — p)/ (2-7r) 3 / 2 6(l + p) for the majority spin species 
and 1 + a[2 - (^ F |) 1 / 2 ]/(27r) 3 / 2 & for the minority spin 
species. This reduction in the effective chemical poten- 
tial and corresponding fall in local density can be under- 
stood in terms of an increase in the local pressure within 
the cloud due to the repulsive interactions between the 
atoms. This pressure inflates the cloud causing the RMS 
radius to rise through a factor of 1 + a(l — p)(3yT+~p — 
yT— ^)/2 5 / 2 tt 3 / 2 &[(1 + p j3/2 + (1 _ p )3/2]. Having ana- 
lyzed the weakly interacting regime it is natural to also 
examine the strongly interacting limit. Here the atomic 
gas is fully polarized so = and n+{r) = 7+ — V(r), 
meaning that the system is firmly in the Thomas Fermi 
regime. We again require that the number of particles is 
conserved which sets the majority spin Lagrange multi- 
plier to rescale by a factor of 2 1 / 2 from its original value if 
there were no population imbalance. Consequentially the 
enhanced Fermi degeneracy pressure dilates the RMS ra- 
dius of a cloud with zero population imbalance by a factor 
of 2 1 / 4 . The key limits of weak and strong interactions 
hold true whatever the true theory of ferromagnetism so 
provide two valuable handles for potential experiments. 

B. Exact analysis of trapped behavior 

Having completed the overview of the trapped behav- 
ior we now turn to consider the ramifications of fluctu- 
ation corrections and self-consistcntly solve Eq. ([5]) for 
the chemical potentials /i±. We then integrate over the 
trap to extract the full behavior of the experimental ob- 
servables, namely cloud size, kinetic energy, and three- 
body loss rate, which for the mean-field limit are shown 
in Fig. [3] The same calculation repeated for fluctuation 
corrections is shown in Fig. [4] A useful reference through- 
out will be the complementary analysis in three dimen- 
sions El- The orthodox Stoner mean- field theory 
predicts that at the onset of ferromagnetic ordering the 
system immediately fully polarizes across the entire trap 
at a/b = V27T 3 / 2 w 7.874, whereas fluctuation corrections 
allow the cloud to adopt partial polarization over the win- 
dow of scattering lengths 3.953 a/b 4.048. In two- 
dimensions as the entire gas polarizes at the same interac- 
tion strength striking features emerge at these respective 
interaction strengths. Current experiments Q can probe 
scattering lengths to ~ 10% accuracy, therefore in cur- 
rent experiments the fluctuation corrected transition will 
also appear to immediately give complete polarization. 

To develop our intuition we first examine the projected 



Mean-field 




FIG. 3: (Color online) The variation of (a) cloud size, (b) ki- 
netic energy, and (c) atom loss rate on ferromagnetic ordering 
with increasing scattering length a/b for the orthodox Stoner 
mean-field theory case. The thin blue dashed line highlights 
the small a/b behavior. The solid lines are at zero popula- 
tion imbalance whereas the dotted line is with an imposed 
population imbalance of 0.5. 

cloud size. In the mean-field limit (Fig. |3Ja)) with weak 
interactions the RMS radius grows linearly with scat- 
tering length as the atoms repel each other within the 
trap. The radius grows following the universal scaling de- 
scribed above. Population imbalance causes the cloud to 
have an initially larger radius due to the increased Fermi 
degeneracy pressure. The entire cloud becomes fully po- 
larized at the same scattering length, a/b = \/2n 3 / 2 , and 
at this point the cloud size immediately adopts its fi- 
nal inflated radius R/Rq = 2 1 / 4 , maintained by Fermi 
degeneracy pressure. This is in contrast to the three- 
dimensional case [ill d2 in which the transition takes 
place over a range of interaction strengths, thus making 
the transition less distinct. Fig. |U[a) shows that fluctua- 
tion corrections drive the cloud expansion faster, causing 
it to dilate rapidly. In contrast to the three-dimensional 
case [l(| HH, this pressure cannot drive the cloud to 
grow larger than the fully polarized size 2 1 / 4 i?g MS . As 
the interaction strength is unaffected by the density of 
atoms, the transition occurs at the same scattering length 
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Fluctuation corrected 




p = o 

p = 0.5 
Qall a/b 



■ (b) Kinetic energy 
0.3 - 




FIG. 4: (Color online) The variation of (a) cloud size, (b) 
kinetic energy, and (c) atom loss rate on ferromagnetic or- 
dering with increasing scattering length a/b when fluctuation 
corrections are taken into account. The thin blue dashed line 
highlights the small a/b behavior. The solid lines are at zero 
population imbalance whereas the dotted line is with an im- 
posed population imbalance of 0.5. 



a w 3.9546 seen in the uniform case. 

The total kinetic energy is probed experimentally by 
releasing the atoms from the trap and imaging them fol- 
lowing a ballistic expansion. Starting from the mean- 
field analysis in Fig. [3Jb) , at zero interactions an initial 
population imbalance increases the kinetic energy due 
to the enlarged majority spin Fermi surface by a fac- 
tor of [I + (j^|) 3 / 2 ]/2. The weak interactions dilate 
the cloud, causing local density and kinetic energy to 
fall with the universal scaling 1 — a(l — p)(3y r T+~p + 

y T -T^)/ 2 5/2 7r 3/2 6 [ (1 + p j3/2 + (1 _ p)3/2]_ Whcn thc 

scattering length is increased beyond a/b = \Z2n 3 / 2 thc 
entire gas becomes ferromagnetic and thc atoms all enter 
the same Fermi surface. This Fermi surface is inflated 
and the kinetic energy plateaus at the final value that is 
2 1 / 2 times that for the non-interacting gas. When fluctu- 
ation corrections are taken into account one recovers the 
variation of kinetic energy shown in Fig. EJb). The fluc- 
tuations drive the transition to take place at a reduced 



interaction strength a w 3.9546 seen in the uniform case. 

The atom loss rate due to three-body recombination 
is T = T (a/b) 6 J n+(r)n_(r)[n + (r) + n_(r)]d 2 r [p. In 
the recent experiment [9( the three-body loss was signifi- 
cant and forced the experiment to be performed rapidly 
and out of equilibrium, and here we study the situation 
in two-dimensions. We start by examining the mean- 
field limit in Fig. [3Jc) , which shows the three-body loss 
integrated over the entire trap. At weak interaction 
strengths the loss rate rises rapidly as T = Fo(a/6) 6 /j, 4 (l — 
p 2 )/87r 2 w. At a scattering length a/b = \[2tv'I' 1 the gas 
across the entire trap becomes fully polarized so n_ = 
and therefore the three-body loss is completely cut off. 
This immediate elimination of loss contrasts the three- 
dimensional case where loss remains until high interac- 
tion strengths, where it forces the experiment out of equi- 
librium [10| | , and also renormalizes the effective interac- 
tion strength Fig. [3] highlights how these effects are 
reduced in the two-dimensional case which could aid with 
the positive identification of the ferromagnetic phase. It 
can also be seen that population imbalance reduces atom 
loss primarily through reduction of the n + (r)n_ (r) term. 
Having studied the mean-field limit we now look at the 
impact of fluctuation corrections on three-body loss in 
Fig. life). The fluctuation corrections drive the ferro- 
magnetic transition to take place at a reduced scattering 
length of a/b = 3.954. This in turn means that the peak 
three-body loss (cx a^ rit ) is significantly reduced. This fall 
in loss rate will mean that an experiment searching for 
signatures of ferromagnetism can be performed nearer to 
thc equilibrium regime which should yield clearer results. 



V. DISCUSSION 

In conclusion, on the repulsive side of the Feshbach 
resonance coupling of transverse magnetic fluctuations 
drives ferromagnetic ordering first order. We studied the 
specific variation of three experimental signatures of fer- 
romagnetism: cloud size, release energy, and atom loss 
rate. Thc formalism highlighted the benefits of studying 
ferromagnetism in two rather than three dimensions. In 
two-dimensions thc effective interaction strength is inde- 
pendent of density and therefore radius in the harmonic 
well. As the interaction strength is ramped upwards the 
entire gas will enter into the ferromagnetic phase at the 
same Feshbach field, whereas in three-dimensions the gas 
first enters the ferromagnetic state at the center. There- 
fore the signatures of the ferromagnetic phase are en- 
hanced in two-dimensions, which should aid the exact 
characterization of the state. At weak interactions these 
observables displayed universal scaling, and the variation 
with an imposed population imbalance was also consid- 
ered. 

One intriguing possibility opened up by the new for- 
malism developed to study fluctuation corrections is fer- 
romagnetic reconstruction into a spin textured state, in a 
matter analogous to the FFLO state in superconductors. 
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This has already been shown to be possible in three di- 
mensions Q and, with enhanced Fermi surface nesting in 
two dimensions, poses an interesting direction for future 
research. 
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FIG. 5: The re-parameterization of the momenta ki,2,3,4. 612 
represents the angle between ki and k 2 . The two momenta 
P12 = ki +k2 and P34 = k3 +k 4 are constrained to be equal, 
P12 = P34, by the Dirac delta function in Eq. (|A1[) . 



Appendix A: Computing the momentum space 
integral 

An important integral Eq. (U) encountered in this pa- 
per has the form 



dk 1 dk 2 dk 3 dk 4 F(fc 1 , fc 2 , fc 3 , fc4)5(k 1 +k 2 -k 3 -k4 



(Al) 

To evaluate this integral one could substitute k 4 = 
ki +k2 — k3, and then integrate over the three parameters 
representing the lengths of vectors ki, k 2! and k 3 , and 
a minimum of three relative angles between these vec- 
tors, giving a total of six integration parameters. How- 
ever, since numerical integration generally becomes pro- 
hibitive with increasing number of dimensions we outline 
a scheme that takes advantage of the fact that the func- 
tion F depends only on the magnitude of the momenta to 
perform the angular integrals and leave a numerical inte- 
gral over just the four vector lengths. A similar scheme 
has been developed in the three-dimensional case 0j- 

The integral is re-parameterized according to Fig. [5] 
The angular integral associated with vectors ki and k 2 
is J w 27rfcifc 2 d#i 2 , where #12 is the angle between ki and 



k2. We now change the variable of the angular integral 
over 9 12 to the vector p 12 = k x +k 2 through the relation- 



(kf 



k\ 



p\ 2 )/2kik2 and so d0\ 



ship cos 9i2 

dp 12 8rk 1 k m2 [4kik 2 2 -(kl+kl-pl 2 f}- 1/2 - This 
expression, and an analogous one in P34 = k3+k4, allows 
us to rewrite the original integral Eq. (|A1[) in terms of 
the parameters p 12 and P34. The momentum conserva- 
tion requirement is imposed by <5(k!-|-k 2 — k 3 — k 4 ) which 
now introduces a new conservation law 8(jp± 2 — P34). This 
sets the two integration parameters equal, p 12 = P34 = p, 
so there is just one integral over parameter p remaining, 
and since the delta function constrains the angle between 
Pi 2 and P34 we must also divide by the phase space as- 
sociated with the angular integration of 2irp. We then 
obtain 



frrf /■min(ki+k2,k 3 +k i ) 

32?r / / / / dfcidfc 2 dfc 3 dfc 4 / dp x 

JJJJ Jmax(\k 1 -k 2 \,\k 3 -k i \) 

F(ki, fc 2 , fc 3 , fc4)fcifc 2 fc 3 fc4p 

^/Ak'fk'i - (fc 2 + kl- p 2 )V4fc^| - + k\- p 2 ) 2 



Finally, we note that the integral over variable p is Carl- 
son's standard elliptic integral of the first kind, which we 
denote by Rp. This yields the final result 



32tt 



d/cidfc 2 d/c 3 d£;4-F(£;i, fc 2 , fc3, ^4) x 



e(fci + fc 2 - 


- |*3 - fc4|)e(fe + 


ki — \ki 


-k 2 \) x 




Rf (0, 1 + 


[(fe 1 +fe 2 ) 2 -(fc 3 -fc4) 2 


[(fci-fe 2 ) 2 


~{k 3 +k 4 ) 2 ] 


•0 


[(fci+fcaP-tfca+fci) 2 


pi-fc 2 ) 2 


-(fc 3 -fc4) 2 ] 



VIKfci - k 2 ) 2 ^ Wa - fc 4 ) 2 ][(fci + k 2 y - (k 3 + kif 

The term introduced to compensate for the angular in- 
tegrals can be efficiently computed by a suitable nu- 
merical library. This four-dimensional integral is now 
better suited to computational evaluation than the six- 
dimensional form of the original expression Eq. (|Aip . 



